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Abstract. This article reviews the main features of the transport properties of 
single fractures. A particular attention paid to fractures in geological materials which 
often display a roughness covering a broad range of length scales. Because of the 
small distance separating the fracture walls, the surface roughness is a key parameter 
influencing the structure of the void space. Studies devoted to the characterization of 
the surface roughness are presented as well as works aimed at characterizing the void 
space geometry. The correlation of the free space is found to be crucially function 
of the failure mechanism (brittle, quasi brittle or plastic.) but also of possible 
shear displacements during the failure. The influence of the surface roughness on 
the mechanical behavior of fractures under a normal load and a shear stress is also 
described. Finally, experimental, numerical and theoretical works devoted to the 
study of the influence of the fracture void geometry on the permeability and on the 
hydrodynamic dispersion of a dissolved species are discussed. 
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1. Introduction 

Since fractures are difficult to avoid and may extend in size over scales ranging from 
microns to the many meters, their flow and transport properties have received a lot of 
attention. Most of the past effort were focussed on fractured in geological materials 
because they are of interest in a number of context ranging from oil or water reser- 
voir management, geotechnical applications such as underground storage of wastes and 
deeper Earth systems such as earthquakes [lj. Amongst others applications, we can cite 
the long term durability of concrete structures or sealing problems in mechanical indus- 
try and more particularly in the design and safety of nuclear power plants or cryotechnic 
rocket engines, where the sealing of some units is performed by a direct contact between 
metallic rough surfaces [21 [3]. 

All the past studies revealed that transport of passive species and mechanical properties 
of single fractures depend almost exclusively on the geometry of their void space pQ. 
Moreover, because existing fractures provide planes of weakness on which further defor- 
mation can more readily occur, fractures also largely control the mechanical behavior 
of materials. Finally, fractures also cause the behavior of large volumes of material to 
be different from that of small scale samples tested in laboratory. 
Past studies clearly point out two important characteristics of flow through fractures: 
The first is that even under normal loads as large as lOOMPa for rocks [I] or QOOMPa 
for pressed metallic joint [2], the permeability of the fractures is still higher than the 
matrix porosity. This demonstrates that the walls of these fracture do not perfectly 
match so that a sufficient free space inside which flow takes place is left. The second 
key feature is the strong channelization of the flow along few preferential paths of low 
hydraulic resistance [HIE]. Even if research efforts were devoted to this subject, the re- 
lation between the void space geometrical characteristics and the fracture permeability 
is still an open question, largely because of the localization of the flow; the latter has 
also a strong impact on the transport of particles or contaminants. In the latter case, 
early breakthroughs of contaminants are often reported for flow in fractured materials: 
the accurate description of such phenomena is still not achieved. 

In this paper, we report briefly recent developments together with the most important 
past results concerning the description of the transport properties of single fractures. 

The first section of this paper deals with the statistical properties of the aperture 
field of single fractures either made artificially or sampled on field sites. Although the 
aperture field have a very variable topology, they can be characterized by a correlation 
length. Its value is found to be crucially related to the genesis of the fracture: For 
instance, plastic deformation and small scale disaggregation following failure prevent 
the fracture walls from mating perfectly at small length scales, in this case fractures are 
characterized by a small - compared to the fracture size - correlation length reflecting 
past physical processes. Small displacements between the fracture walls also prevent 
the fracture from mating. We will show that because of the surface roughness that 
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may extend over scales comparable to the fracture size, the aperture is anisotropic 
with a large correlation length in the direction perpendicular to the shear displacement 
direction. The evolution of the void space as well as the geometry of the contact zones 
are reported for various stress conditions (normal load or shear stress). 
The void space is a key parameter controlling the transport property of fractures. For 
instance, the amount of water that may flow under a pressure gradient through a fracture 
is crucial in many applications. Yet, all past experimental observations found that this 
quantity cannot be simply related to the geometrical parameters such as the mean 
fracture aperture. This problem as attract a lot of attention but no satisfactory model 
as been yet found. 

The heterogeneous nature of the aperture also strongly influences the transport of solute 
particles: because of preferential channels inside which a large part of the flow takes 
place, an important fraction of solute is transported by convection over large distances 
rapidly while the rest of the solute lags behind resides in reduced velocity regions. As a 
result, classical description of dispersion using a Fickian approach breaks down at the 
fracture scale. The studies performed on that topic are described in the last section 
(Sec. |5]) of this paper. 

A key parameter controlling the aperture is the roughness of the fracture walls. The 
following section discusses results concerning the statistical description of the roughness 
of fracture surfaces. Past results on single fractures propagating in disordered systems 
such as rocks report that the roughness increases non linearly with the fracture size. 
We will then show how this multiscale roughness influences the transport properties of 
fractures. 

2. Roughness of the fracture walls 

The structure of the aperture field is strongly related to the relative position of the 
fracture together with the statistical distribution of the asperities. Naturally occurring 
rough surfaces are characterized by a roughness present over a broad range of length 
scales. As illustrated by figure HJ these correlations typically manifest themselves 
through the increase of the amplitude of the roughness with the observation scale 
considered. Statistical analysis of the scaling dependance of the surface roughness on 
the fracture size reveals that the roughness of fractured surfaces is well described as 
a self-affine fractal structure 0. This means that the spatial distribution of the local 
heights h(r) is statistically invariant under the transformation r — > Ar and h — > X^h 
where (, the self-affine or Hurst exponent, ranges from to 1. So, giving the correlation 
function 7(5) = (Ah 2 ) to find a height difference Ah = h(r) — h(r — S) over an horizontal 
vector S, the scaling invariance leads to: 

7 (A<5) = \ 2 h(S). (1) 

| Numerous works have characterized the surface roughness through fractal dimensions D. The values 
reported agreed with the inferred values of £'s; henceforth no direct reference to D will be made. The 
difference between this two exponents is discussed in [7]. 
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Figure 1. Map of a fractured surface of a granite block broken under a tensile stress. 
The heights h(x, y) were recorded by a mechanical profilometer surface. The size in 
the [x, y) plane of the surface is 200 mm x200 mm, over this area rms of the heights 
is 3.5 mm and the height difference between the highest and lowest points is 20 mm. 

Figure [2] shows the variation of the function j-fl computed for distances S oriented along 
(circles in Fig. [2]) and normal (squares in Fig. [2]) to the direction of propagation of 
the fracture front. For both direction, 7 is a monotonously increasing function of the 
lag distance \5\ considered indicating a clear increase of the surface roughness with its 
extension. In addition, the increase is linear indicating that 7 follows a power law of 
5 with an exponent ( = 0.75 ± 0.05 which shows only a weak dependance with the 
orientation considered. Exponents close to 0.8 have been reported in a large number of 
previous studies and the hypothesis of its universality was furthermore conjectured [8j. 
More recent works indicate that this exponent may take two different values, depending 
on the type of fracture materials: 0.8, for most of the materials including glass, cements, 
granite, tuff [HI dU [12j H31 HH US, EH] while ( ~ 0.5 for materials such as sandstone 
[T7] [18] , calcite [19] or materials with similar microstructure such as sintered glass beads 
blocks [20]. The geometry of the microstructure is not the sole parameter controlling 
the roughness. Other experiments performed using sandstone reported a roughness 
exponent close to ( = 0.75 ± 0.03 [21] disagreing with the observation of [T7J ITS] . 
However, in this work [21] a fracture process zone measured by acoustic emission ahead 
of the fracture front was actually present while in the latter experiments fTTJ [18] such 
zone was not existing. The existence of damage zone where brittle micro cracks spread 
affecting the growing of a single brittle crack was then found to be a key mechanism 
explaining the existence of two types of fracture roughness |21j. These non linear and 
plastic effects were further conjectured to correspond to a Family Vicsek self-affme 
morphology |22j. 

The effect of the fracture size on its roughness is a key issue which was studied by few 
authors: Matsuki et al. [16] performed measurement over metric scale granite fractures 
that do not revealed any variation of the exponent ( with the fracture size. The same 
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Figure 2. Variation of the log of the correlation function 7 (mm 2 ) as function of 
the log of the lag distance in \S\ (mm) for vectors oriented along the direction of 
propagation of the fracture font: circles and squares are for measurements normal to 
that direction. The dashed line has a slope of 1.5. The surface studied is the granite 
surface of size 200 x 200 mm displays in Fig. [T] [60 . 

conclusion was drawn from the field measurement performed by Schmittbuhl and co 
workers [13]. Yet, by combining millimetric scale observations with measurements at 
the metric scale, Brown and Scholtz [23j still observed that the rms deviations of the 
height of the fracture walls increases with their length but, this time, with an exponent 
decreasing with the size of the fracture. 

3. Dependence of the aperture field on mechanical deformations 

3.1. Aperture field reconstruction 

The topography of a single surface does not tell the full story: the key parameter con- 
trolling mechanical and transport properties of the fracture is the correlation between 
the roughness of the two fracture walls. Different methods have been used to map the 
fracture void space. In the most popular method, the map is reconstructed by using the 
maps of the heights of the two individual fracture walls as measured by profilometry 
0, [2U EH]- This method usually requires to adjust the halves in order to get the best 
match estimates by, for example, minimizing the rms of the aperture. The resulting 
field corresponds to the aperture of the fracture under zero normal stress and with the 
two walls as close as possible to a perfect contact over their full surface. The need of a 
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reconstruction procedure can be avoided by performing injections of fluorescent epoxy 
into the void space [251 EH]. In this case, the aperture is reconstructed through image 
analysis of sections; for practical reasons the section are large making difficult the frac- 
ture reconstruction. Measurements were also done an transparent replicas of natural 
fractures [27J ES EH]. Dyed fluid is then injected into the fracture void and the local 
light absorption is related through an appropriate calibration law to the local fracture 
aperture. Recent development of X-Ray tomography and the availability of very intense 
beams offers a new and efficient way to measure the void structure of porous and frac- 
tured media with a spatial resolution that can be as good as 5/xm on centimeter scale 
objects [IHl EH1 EH [32] . Finally, technical improvements together with the development 
of new surface profilometry techniques (such as atomic force, scanning electron or con- 
focal microscopes and laser scanner) allows now the acquisition of 3d maps of the halves 
of the fracture [U [161 EE1 [33] with a very good precision. 

With most of these techniques it is possible to determine the statistical distributions 
of the apertures as well as their spatial correlations. These distributions are usually 
well fitted by normal [HI [2H [261 EI] or log-normal [281 El] distributions with mean of 
the order of few hundred microns and standard deviation of the order of 0.1/zm under 
zero normal load. Distributions adjusted by log-normal distribution law are usually 
observed on fractures that have been submitted to a normal stress: in this case, the 
aperture distributions are shifted towards lower values as the surfaces come in close 
contact: normal distribution may then with increasing stress and contact be fitted by a 
log- normal distribution [H ESI El]. The characterization of the correlation through the 
adjustment of the semivariogram 7 of the aperture field reveals that most of the field are 
correlated over length scale ranging from few hundred of microns up to few millimeter 
independently of the rock type 0HIS1I231ISIII11I2ZIEIII2S1E31I2HII23- 

3.2. Mechanical deformation and aperture variations under normal load 

When a normal stress is applied to a preexisting fracture, its aperture reduced and 
contact zones appear resulting in a non linear variation of the deformation with stress: 
the rate of deformation is highest at low stress values which reflects an increase of the 
fracture stiffness. After repeated cycling these non linearities which reflects changes in 
the void space geometry are still observed [351 E3]: some voids get closed and addi- 
tional surfaces come in contact. The details of the void space geometry were studied by 
Pyrak-Nolte and co-workers [35] by means of an injection of a low melting-point alloy 
allowing one to obtain a cast of the void space. Pictures of the two surfaces are then 
taken and digitized. The two maps obtained in this way show the locations where metal 
has adhered on each surface, but do not give the value of its thickness. The geometry of 
the area in contact was also studied by Gentier [36] who placed thin film of plastic sheet 
between the fracture walls. When squeezed, the optical quality of the film is locally 
altered: image analyzing is then used to determine the region where the fracture walls 
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are in contact. Experiments repeated under normal stress up to 80 MP a indicate that 
contact areas are heterogeneously spread over and represent a total area of the order of 
60% of the fracture area jU [361 EH]- Recently Sharifzadeh et [33] report for a fractured 
granite block contact area fractions of the total surface corresponding to the contact 
area ranging from 86% for a normal stress of 1 MPa up to 99% at 10 MPa. The latter 
value indicated a near perfect match between the two surfaces. 

3.3. Mechanical deformation under shear 

Recent experimental investigations were intended to determine, visualize and interpret 
aperture distribution under various shear displacements and different normal stress vari- 
ation histories. All these works demonstrate that the application of a shear stress alter 
significantly the void space. Two main features were observed: a dilatancy of the frac- 
ture resulting in an increase of its aperture and in the appearance of damage zones. 
The relative magnitude of the two phenomena is controlled by the normal stress. If the 
latter is low compared to the rock strength than gouge material is not be produced and 
protuberant asperities behave like rigid material [33J. The opposite situation - normal 
stress overcoming the strength of the material - is achieved by, for instance, using ce- 
ment mortar cast of rock fractures [371 EE]. 

When a shear stress is applied to one of the fracture walls, small variations of the fracture 
aperture are first observed up to a shear displacement of a few millimeter [10, EHJ [371 133] 
because of the motion is stopped when the asperities of the upper and lower surfaces 
come in contact. This is followed by an uplift of the two halves reflecting the dila- 
tion, as the shear stress increases simultaneously up to a peak. After that, the shear 
stress rapidly decreases to a residual value and becomes almost constant. At a critical 
point near the peak value of the shear stress, for which the latter becomes equal to 
the strength of the asperities in contact, all the asperities in contact are simultaneously 
sheared. As a result, the aperture distribution suddenly changes because the roughness 
of the two surfaces do not match any more. Increasing further the shear results in a 
further increase of the aperture and in a reduced contact area. For high strength mate- 
rials, the fraction of the area where the surfaces are in contact drops abruptly when the 
peak value of the stress is reached [331 H]- Meanwhile, the mean aperture increases by 
more than an order of magnitude [33J corresponding to an increase of the permeability 
by three orders of magnitude [4]. For more brittle materials, damage is initiated in 
areas of the fracture with the steepest slope and their extensions in the fracture plan 
are strongly related to the amount of vertical displacement [37J. The gouge particles 
created during the shear displacement may in turn influence the fractures properties [38J . 
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3.4- Modifications of the aperture field under shear 

When the fracture halves are brought together so they are separated by a minimum free 
space, the aperture frequency distribution follows a normal or a log-normal distributions 
and the regions in contact are evenly distributed. When the shear increases, the aperture 
distribution is shifted towards the high aperture values indicating an increase of the mean 
fracture aperture. In the same time, the distribution broads and becomes well adjusted 
by a normal law [H HO]. Under such conditions, the contact areas are then spatially 
localized and an heterogeneous distribution of the aperture is observed. Moreover, 
the analyze of the spatial correlation shows that the aperture field is characterized by 
a strong anisotropy with a correlation length which is much larger in the direction 
perpendicular to the displacement [U [27] . 

This crucial shift in the aperture field is illustrated quantitatively by considering 
the simplified configuration in which fractures made of complementary walls that are 
separated both vertically in order to open the fracture and laterally in order to mimic 
shear displacements (if any). Under such conditions, if a is the mean aperture and u the 
relative shear displacement in the plane of the fracture, the local aperture at a point f 
satisfies : 

a(r) = a(r) — a(f + u) + a. (2) 

This equation allows one to determine the aperture field a(r) for a given shift u from 
the experimental surface profile map a(r) assuming that the surface match perfectly 
for zero normal and lateral shifts. The Fig. [3] shows examples of aperture fields com- 
puted from the heights map of a granite fracture. The size of the fracture is 60 mm 
and lateral shifts representing only 4% of the fracture size were applied along the two 
diagonal of the fractures. Even if the displacement is small, the aperture fields show a 
remarkable anisotropy with a strong correlation in the direction perpendicular to the 
displacement. Such an anisotropy was also observed by Matsuki and co-workers [16] in 
the case of two walls which do not match perfectly under zero shear displacement as 
soon as a shear displacement large enough to generate heterogeneities with a correla- 
tion length overcoming the correlation length observed in the absence of displacement 
was imposed. The correlations in the aperture field can be characterized through the 
correlation function j(8), defined in Sec. [21 which measures the spatial correlation of 
the aperture field between two points separated by the 2D vector 8. The orientations 
of the lag 8 respectively parallel and perpendicular to the shear displacement u are of 
special interest (in the following, the lags are respectively referred to as || and _L). Fig.H] 
displays the semivariances 7(5ii) and 7(<5j_) of the aperture fields display in Fig. [31 For 
small magnitudes 8 and for a given exponent, the apertures are very similar (high spa- 
tial correlation at short distances) and 7 is small. As \8\ increases, so does 7, because 
the spatial correlation of the aperture decreases. If \8\ increases further, up to a value 
higher than the correlation length of the aperture field, the apertures become uncor- 
rected and 7 should reach a same saturation value for both orientations of 8 equal to 
twice the aperture variance <j 1 a =< (a(x, y) — a) 2 > x ,y This value corresponds to a ratio 
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Figure 3. Aperture field obtained by displacing a surface height map of a granite 
fracture according to Eq. [2] The magnitude of the shear displacement is in both cases 
u/L = 0.04 where L = 60 mm is the fracture size. The shear directions are indicated 
by arrows. Figures (c) and (d) represent the aperture for a threshold value set to half 
the mean fracture aperture (white: a < a/2, black: a > a/2). 

7 /2o~1 = 1 (dashed line in Fig. 0j) which is actually the asymptotic value observed at 
high \5\ values on Fig. HI Even though the asymptotic values are the same, the way 
in which they reach this saturation value differs greatly depending on the orientation 
of S with respect to the shift. The normalized semivariances corresponding to the per- 
pendicular direction, never exceed the saturation value and reach it in an overdamped 
way. On the other hand, the functions 7/2<r^ for \S\ parallel to the shear displacement 
display a small overshoot, reflecting a local anticorrelation. These different behaviors 
result directly from the large-scale anisotropic structures observed in Fig. [3l 

The history of the stress variations experienced by the fracture has a strong influence 
on its aperture. Under high normal load, the contact area between the two halves is 
important and the contact zones are distributed heterogeneously in the fracture plane. 
As a consequence the flow and the geometry of the path of least resistance will be greatly 
influenced by the fracture spatial distribution of the void space. When the fracture walls 
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Figure 4. Variation of 7 normalized by 2a\ where a a is rms of the aperture as 
function of the lag \6\ normalized by the value of the shear displacement u for direction 
perpendicular (filled symbols) and along (open symbols) the shear displacement. The 
fields are shown in Fig. [3J circles and squares correspond respectively to field (a) and 
(b) 



have undergone, in addition, a shear displacement, high aperture channels develop in the 
direction normal to the shear and represent paths of reduced hydraulic resistance. The 
following section reviews some of the studies and the most important findings concerning 
the flow of a single liquid in a fracture. 

4. Experimental and numerical estimation of the fracture transmissivity 

As a first approximation, it is natural to view fractures as made of parallel flat plates 
with a separation equal to the average mean distance separating the fracture walls. In 
this case, the macroscopic transmissivity is proportional to the cube of this aperture. 
Yet, a large number of experimental measurements of the transmissivity of fractured 
rocks reveals that the transmissivity decreases more rapidly than the cube of the mean 
aperture [TJ. These findings demonstrate that as soon as the fractured surfaces are in 
contact (or close to contact) the complex geometry of the structure of the free space 
where flow takes place strongly alter the flow, leading to an increase of the hydraulic 
resistance. This non linear relation between the cube of the fracture aperture and the 
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transmissivity results from the strong localization of the flow along a few preferential 
paths which is found to increase with the degree of heterogeneity. A classical way to 
study the coupling between the fracture geometry and the flow structure is through the 
Reynolds approximation which - by neglecting the out of plane flow movement - reduces 
the 3D flow problem into a 2D problem. The next section summarizes the hypothesis 
used to derive the Reynolds equation from the Navier-Stokes in the context of flow in 
fractures. 



4-1. Equation of motion of the fluid motions and definition of the fracture 
transmissivity 

Three elements are necessary to describe the motion of a fluid: the conservation of mass 
and momentum and the boundary conditions. For an incompressible fluid (this is the 
case for most of the liquids), the mass conservation equation reduces to V.v = where 
v is the velocity vector describing the motion of the fluid particles. For a Newtonian 
fluid, the momentum equation reduces to the Navier-Stokes equations given by [41 J: 



p— + (v.V)v = -Vp + pS7 2 v 1 (3) 

where p, p and p are respectively the dynamic viscosity of the fluid, its density and the 
total pressure. 

In most flows through fractured media, the effect of inertia can be considered as 
negligible to that of the viscosity. This can be made quantitative through the estimation 
of the Reynolds number Re given by: 

Re = ^ (4) 

Lit 

where I is a characteristic distance over which the aperture varies noticeably in the 
direction of the flow [48]. When the inertia term is neglected (Re << 1) and in the 
steady state regime, the Navier-Stokes equations can be reduced to the so called Stokes 
equations [41] : 

- Vp + pN 2 v = o. (5) 

This assumption breaks down in locations where the apertures displays fast variations: 
then, the effect of the inertial must be taken into account and, in regions where the 
velocity decreases abruptly, recirculation eddies appear. These effects result in an 
additional dissipation which influences in turn the velocity field, leading to the reduced 
permeability observed experimentally (6j E2J [63], 154] . However, recent experimental 
observations [65] demonstrated that such non linear flow regimes are observed at the 
fracture scale only for Reynolds number greater 10. 

Together with the mass conservation equation, the three components of the Stokes 
equation form a set of 4 coupled equations and that can be solved analytically only for 
very simple cases. The most widely used hypothesis reducing this 3D problem to a 2D 
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one is to assume that the tortuosity of the flow lines across the local aperture (i.e. along 
z in Fig{T]) is negligible. This is equivalent to consider that the pressure variation across 
the gap is small in regard to the pressure gradient in the fracture plane i.e. §f << f 2 
(lubrication approximation). As a result, only the components of the velocity vector 
aligned with the pressure gradient are non zero. Together with a dimensional analysis 
of the variation of the velocity with respect to the direction |48j, one show that Eq.© 
reduces then to: 

Vp(x,y) = fx—. (6) 

Integrating this equation with respect to z bearing in mind the no slip boundary 
conditions at the fracture walls yields to a velocity profile which is parabolic (the velocity 
being null at the walls and maximum half way between them). 

The integration of the velocity profile in the fracture gap leads to a estimation of the 
flow rates per unit of fracture width Q x and Q y which are given by: 

Q x = -T(x,y)^ (7) 

Q y = -T(x,y)^ (8) 

where T(x,y) = * s the local fracture transmissivity in m 2 .s~ l . Together with the 

mass conservation, V.Q, the previous set of equations reduce to the Reynolds equation: 

V.(T(x,y)Vp) = (9) 

where the derivatives are taken over the (x, y) plane. 



4-2. Transmissivity of natural fracture and its modeling 

Quantitative comparisons between experimental measurements and numerical estima- 
tions were rarely performed in configurations in which the detailed of the geometry of 
the void space was known. Most of these works obtained the macroscopic transmissiv- 
ity by solving the Reynolds approximation (See Eq.^SJ). This approximation, even if it 
usually leads to an overestimation of the flow velocity or equivalently of the hydraulic 
aperture that increases when the fracture aperture is reduced (23, [29j [42], [43], is used 
extensively to deal with flow problems in fractures |14| . Mourzenko et al [H] reported 
systematic comparisons between numerical solutions for flow in rough fractures obtained 
using the Reynolds approximation (see Eq.flH}) and the 3D Stokes equations given by 
Eq.Q. The simulations are found to give similar results only for wide and smooth 
fractures; but, when the fracture aperture is decreased or the aperture fluctuation is 
increased important differences between the two sets of results are observed. Compa- 
rable conclusions were obtained by Brown et al [15] when comparing the two types of 
solutions for flow in sinusoidal channels. The geometrical validity of the application of 
the Reynolds approximations is discussed in |34j. 
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Insofar flow through the fracture is modeled in the framework of the the lubrication 
approximation, the problem of finding the macroscopic transmissivity of the fracture 
reduces to the problem of finding the effective conductivity of an heterogenous 2D 
conductivity field. Many theories have been developed in the past years to estimate this 
quantity. Among the most popular is the small perturbation analysis [46J. In the context 
of fracture flow, this approach allows one to estimate the effective transmissivity for 
any aperture distribution - including normal, lognormal distribution [37], unidirectional 
ripples jlH] - by means of a first order expansion given by [3Tj: 



where 9{u) is the Fourier transform of the correlation function 1 — , ~f(5)/2a^ defined in 
Sec. [31 Equation ffTOl clearly demonstrates that while the magnitude of the macroscopic 
transmissivity is controlled by the large aperture zones that control a, the less open 
zones also play an important part by reducing the transmissivity of the fracture when 
it is closed. 

The prediction of Eq. tflOl) compare reasonably well with several numerical simulations 
[501 [10] an d experimental results [H]. Yet, this perturbation theory ignores the possible 
channelization of the flow [HI G2E [H7J [HE] , so that deviation from this theory with exper- 
imental observation are expected when the heterogeneity of the aperture is enhanced, 
for instance, when large normal load are applied. 

Fluid flow through fracture is also hindered by the presence of contact areas which, 
under a large normal load represent an important fraction of the total fracture area. For 
fraction lower than 25% and under the assumption that the fracture aperture is constant 
outside the contact zone so that the governing equation for fluid flow is a Laplace equa- 
tion, Zimmermann and co workers [51] obtained an analytical expression of the effective 
transmissivity for various simplified shapes of the contact area (circular or elliptical). 
For contact area occupying a large fraction of the surface, flow is controlled to a large 
extent by the location and distribution of the constrictions connecting the larger voids 
into continuous paths [52J. It is them tempting to map the void structure of tightly 
mated fractures walls onto a discrete network. Recently Plouraboue et al [3] developed 
an efficient method based on mapping the continuous transport equation onto a discrete 
network of conductance function of the local geometry of the fracture aperture [53J. 

4-3. Transmissivity variations induced by a shear displacement 

As discussed in Sec. EJ shear displacement causes the fracture to dilate resulting in an 
increase of the fracture transmissivity [40J . The displacement also induces an anisotropy 
of the aperture field. The correlation length of the aperture is indeed much larger in 
the direction normal to the shear displacement then parallel to it. As a result, the 
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transmissivity is also anisotropic with a larger value in the direction perpendicular to the 
shear. More precisely, numerical [U [HI El] and experimental [U [161 CHI [59] investigations 
report an increase of the anisotropy with the shear displacement with a transmissivity 
increases in the direction normal to u while it drops off parallel to u. When the fracture 
aperture is large compare to the aperture fluctuation (i.e. small a a /a), Auradou and 
coworkers [15] found that the behavior is adequately reproduced by considering the 
fracture as made of parallel channel along with the aperture is constant. For this model, 
in the direction where the aperture constant the transmissivity is proportional to the 
arithmetic mean of the cube of the aperture (i.e. T oc< a 3 >) [18] while in the direction 
where the aperture is varies the transmissivity is proportional to the harmonic mean of 
the cube of the aperture (i.e. T oc< a -3 > _1 ) [48J. Shear displacement also lead to 
channelization flow: but compared to the previous situation, the channelization is due 
here to the structural heterogeneity of the aperture rather than to the broad distribution 
of the local hydraulic transmissivities observed when the fracture get closed. 

5. Hydrodynamic dispersion of dissolved species 

Solute transport is of great interest for instance in fractured rocks because of its 
implication in groundwater pollution, C0 2 or nuclear waste sequestration in geological 
formations, and oil recovery but also in corrosion processes where the transport of 
dissolved species partially controls the leaching rate. At the scale of a single fracture, 
dispersion of the dissolved species results from the combined action of the complex 
velocity field (varying both in the gap of the fracture and in its plane) and of mixing 
by molecular diffusion. The latter allows solute particles to move from one streamline 
to another and homogenizes the spatial distribution of the tracers. In the classical 
approach, tracer particles are assumed to perform a random walk superimposed over a 
drift velocity. This latter is the average of the fluid velocity over an appropriate volume 
(the representative elementary volume or REV) while smaller scale variations induce 
tracer spreading. At the REV scale, the average C(x,t) of the tracer concentration 
over a section of the medium normal to the mean displacement satisfies the convection 
diffusion equation (ADE): 

dC(x,t) TT dC(x,t) n <9 2 C(x,t) 

dt ~ U dx + dx 2 [ } 

where D is the longitudinal dispersion coefficient, U the mean velocity of the fluid (par- 
allel to x). The value of D (or equivalently the dispersivity a = D/U) is independent of 
both time and the traveled distance: it is determined by the combined contributions of 
molecular diffusion (characterized by a coefficient D m ) and advection. The relative order 
of magnitude of these two effects is characterized by the Peclet number : Pe = Ua/D m . 
Experiments on cores containing a single fracture [691 E2 El [68] focused on the break- 
through curves measured at the outlet and analyzing the dispersion from the overall 
fracture. The dispersivities determined on the basis of the traditional Fickian advection 
dispersion models given by Eq. ( Till were found to be close to values inferred from the 
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small perturbation theory [151 [77]. The observations suggested that dispersion is con- 
trolled (like in 3D porous media [71]) by spreading due to velocity variations associated 
to the geometry of the void structure. The latter determines the correlation length of 
the velocity field, leading to the so called geometrical dispersion regime. 
Experiments on transparent cell with walls with a roughness characterized by a small 
correlation length (~ 1 mm compared to pluri centimetric fractures) were them devel- 
oped in order to study the evolution of the mixing front with the distance [7U [72j [75]. 
An example of the concentration fields observed in such fractures is shown in FigJSh- 
The influence of small scale heterogeneities located in the mixing zone (central part of 
the figure) is clearly visible in the concentration field. These elongate streaks result from 
the velocity fluctuations created by the obstacles; they have however of small lateral ex- 
tension so that molecular diffusion is sufficient to smear out these lateral concentrations 
gradients. The quantitative analysis of the concentration front as of its evolution with 
time confirms that the dispersion process is characterized by a dispersivity which does 
not evolve with time (expect close to the fracture inlet). Furthermore at low Peclet 
numbers but larger than 1 (so that pure longitudinal molecular diffusion is negligible), 
the dispersivity is found to be constant and close to the value inferred from perturba- 
tion analysis [75]. However, as the Peclet becomes larger a second regime in which the 
dispersivity increases linearly with Pe is observed. This new regime corresponds to the 
Taylor dispersion mechanism [77] arising from the velocity profile in the fracture gap 
(the velocity is maximum half way between the walls and cancels out at the boundaries.) 
that creates concentration gradient within the gap which is reduced through molecular 
diffusion [SI] . 

However, these experiments conduct in transparent models with a short correlation 
length of the fluctuations of the aperture were not able to reproduce the steps observed 
on breakthrough curves [69j due to localization of most of the flow in a few preferential 
channels. Recently experiments on transparent replica of a real fractures have been 
developed J7SJ [291 [79] in order to study the evolution of the tracer concentration as 
function of the traveled distance in a geometry closer to that of natural fractured rocks. 
For sandstone fractures Bauget and Fourar [29J found that the solution derived from 
the ADE equation appears to be unable to model long-time tailing behavior because of 
the occurrence of flow channels extending along the full length of the fracture. On the 
basis of such observations, models considering that dispersion mainly results from the 
velocity contrasts between such macroscopic channels has been developed [69 l 151 180 1 [29] . 
In this model the fracture is represented by a stratified medium and the transport be- 
tween the inlet and the position of observation is characterized through an heterogeneity 
factor defined as the ratio of the standard deviation of the distribution of the local per- 
meability and of the average permeability of the equivalent medium. Several different 
distributions of the permeability assigned to each strata were considered log-normal [69J , 
normal [29J or computed from the aperture field [80J . 

In the study performed by Bauget and Fourar [29J, the heterogeneity factor is a fit- 
ting parameter obtained by adjusting the breakthrough curves derived from the model 
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with the experimental observations. This factor is found to decrease with the distance 
indicating a decrease of the velocity contrast between the flow lines corresponding to 
an increase of interconnections between the flow paths. This study indicate that tracer 
transport in a single fracture over short distance may be considered as resulting from 
poorly connected parallel channels but as the distance increase the number of possible 
interconnection between these isolated channels increases which slowly reduces the in- 
fluence of the localization of the flow. 

The situation were the fracture walls are laterally displaced is analyzed in Refs [751 179] . 
The corresponding model fractures are characterized by a strong structural heterogeneity 
with channels of least hydraulic resistance developing in the direction perpendicular to 
the shear displacement. When flow occurs along that particular direction, solute parti- 
cles rapidly spread (see the long structures in FigHb) along these channels (white regions 
in Figj5b). A systematic comparison between the geometry of the iso-concentration front 
c(x, y, t) = 0.5 and the stratified model [751 EH] demonstrate, for this flow configuration, 
that the heterogeneity factor constant along the full length of the fracture remains over 
the entire fracture length. These results are not valid any more when the flow is parallel 
to the direction of the shear displacement (i.e. perpendicular to the channels) f79|: in 
this case the breakthrough curves are fitted with the solution of the ADE equations and 
mainly reflect the dispersion induced by the velocity profile in the fracture gap (Taylor 
dispersion [51]). 

In summary, laboratory scale experiments, as well as theoretical investigations, 
demonstrate that the channelization of the flow has a large influence on transport 
processes. However, the importance of this phenomena is related the structure of the 
aperture - a fracture characterized by a broad distribution of aperture will more likely 
display a strongly channelized structure - but also to the flow direction in the context 
of a fracture characterized by a structural heterogeneity. The channelization creates 
concentration gradients transverse to the flow, molecular diffusion tends to reduce such 
inhomogeneities. Thus, if the distance between the channels and/or the Peclet number 
are sufficiently small, transport of solute in the direction transverse to the flow direction 
may be sufficient to reduce the spreading du to the channels: a Fickian dispersion regime 
is then expected to be observed. 

6. Discussion and conclusion 

This paper has discussed recent works dealing with flow and transport in single rough 
fractures. In this context, the statistical properties of the fracture aperture is a key 
parameter. However, fracture aperture is a difficult parameter to quantify and mea- 
sure. The difficulty of these measurements is related to the high precision required in 
order to model the flow field: for instance the permeability of a fracture varies like a 3 
where a is the fracture aperture. Previous investigations have shown that aperture is 
dependent upon stress history, normal displacement, shear displacement and scale of 
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Figure 5. Examples of concentration fields coded in gray level observed on two 
fractures, a) fracture made of a flat plate place above a plate where the roughness 
is obtained by a random distribution of cylindrical obstacles - dots - (See [75] for 
details). The free space above the obstacle is 0.37 mm and the mean fracture aperture 
is 0.65 mm. b) concentration field when flow occurs in a fracture with sheared walls 
along the channelized direction. The arrow indicates the flow direction (See [78] for 
details), c) Aperture field of the fracture with sheared complementary walls. Mean 
aperture: a = 0.75 mm, and shear displacement u — 0.75 mm (see the black arrow). 
The shear is perpendicular to the flow direction. White: regions where a(x, y) overcome 
a/2 (Black is for a < a/2.) 
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study. Under zero shear displacement, the roughness of the two fracture walls may 
not match any more at small scale because of irreversible deformation occurring during 
the failure process or induced by physical interactions with the host fluid. In these 
cases the aperture is randomly distributed and is well described by a normal or a log 
normal distribution. Under normal load, the fracture closes and the areas in contact 
are heterogeneously spread over the fracture plane. During shear, the aperture sharply 
increased at first - fracture dilation - with the shear displacement until a shear peak is 
reached. The aperture distribution then becomes spatially localized. The joint surfaces 
roughness - which is often suitably described by a self-affine multiscale geometry - is 
the main origin of the spatial localization of aperture distribution. For high strength 
materials, the contact ratio decreases rapidly and reaches a small value for large shear 
displacements. For relative brittle materials, asperities of the surface walls are damaged 
during the shear displacement. The damage zones are created in the region where the 
fracture walls have a steep slope creating. The surface roughness is thus - again - a key 
parameter controlling the spatial repartition of these zones. The latter are heteroge- 
neously distributed over the surface. 

The flow of the fluid in a fracture is characterized by a strong localization: the flow 
is concentrated over a small fraction of the fracture area. Channelization arise from two 
mechanisms: the first result from the widely distributed local hydraulic properties while 
the second comes from the structural heterogeneity induces by possible shear displace- 
ment of the fracture walls. The channeling phenomenon has important consequence: for 
instance it imply that solutes transported in channels have potentially far less exposure 
to the rock matrix than in a constant aperture fracture [82J, limiting the reaction, ad- 
sorption, and diffusion of solutes within the rock matrix. In the context of nuclear waste 
management, the ability of the rock matrix to retard solute movement through a frac- 
tured, crystalline rock may be more limited than if all the fracture surface was equally 
available to water flow. A quantitative estimation of the degree of channelization is still 
an open issue, but as point out by past studies, their estimation require information on 
the past hydromechanical history of the fracture. 

Most of the studies reported in the present paper considered small scale 
measurements, a few tens of centimeters in length. Detailed measurements on larger 
scale fractures are mostly limited to measurements of fracture traces [I], they indicate 
the existence of jumps and discontinuities that come, for instance, from changes in 
the loading conditions or merging of fractures. These events create longer length scale 
correlations. Thus at field scales, fracture walls geometry might not appear self-affine 
and the problem of scale effect is still an issue that needs to be addressed in order to 
gain an improved understanding of in-situ transport properties. 
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